Spin polarizability of a trapped superfiuid Fermi gas 
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The polarization produced by the relative displacement of the potentials trapping two spin species 
of a dilute Fermi gas with JVf = iVj. is calculated at unitarity by assuming phase separation between 
the superfiuid and a spin polarized phase at zero temperature. Due to the energy cost associated with 
pair breaking, the dipole magnetic polarizability vanishes in the linear limit and exhibits important 
deviations from the ideal gas behaviour even for displacements of the order of the size of the atomic 
cloud. The magnetic behaviour in the presence of different trapping frequencies for the two spin 
species is also discussed. 
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It is well known that Fermi superfluids cannot be polar- 
ized by an external magnetic field unless the field exceeds 
a critical value. This effect is directly associated with the 
occurrence of a gap in the spin excitation spectrum and is 
the consequence of the existence of pairs. In a dilute gas 
the nature of these pairs depends crucially on the value 
and the sign of the s-wave scattering length. For nega- 
tive and small values of the scattering length the pairs 
coincide with the Cooper pairs of BCS superconductiv- 
ity. If the scattering length is instead positive and small, 
they can be identified with real molecules which, due to 
their bosonic nature, undergo Bose-Einstein condensa- 
tion (BEC). With the aid of Feshbach resonances it is 
now possible to experimentally control the transition be- 
tween the BCS and BEC regimes and to investigate the 
challenging unitary regime where the scattering length 
is much larger than the interparticle distance. The ef- 
fect of spin polarization on these novel quantum phases 
has already been the object of experimental [J, Q, LJ and 
theoretical 0,0.0 investigations. In particular the po- 
larization has been shown to give rise to a phase separa- 
tion between superfiuid and non-superfluid components, 
although the detailed structure of the phase separation 
is still far from being completely understood. 

The purpose of this paper is to investigate directly the 
response of these systems to a spin-dependent external 
field, taking advantage of the fact that the trapping ge- 
ometry suggests a natural way to generate an effective, 
position dependent, "magnetic" field. This is achieved, 
for example, by an adiabatic separation of the external 
potentials confining the two spin species. If the gas is 
nonintcracting the two spin clouds will move rigidly in 
opposite directions giving rise to a spin dipole moment 
per particle 



D(d) 



N 



(1) 



equal to one half of the distance separating the two wells. 
In (JIJ N is the total particle number and d the displace- 
ment of the trapping potential. If instead the system is 
interacting and superfiuid, it will exhibit a resistance to 
spin polarization. In particular the linear response will 
vanish. 

In order to calculate the effects of the local polarization 
induced by the displacement of the confining potentials 
we will make use of the phase diagram of uniform mat- 
ter and of the corresponding equations of state which are 
available in some relevant regimes, the equilibrium be- 
tween the different phases being obtained by imposing 
that the corresponding pressures be equal. In the follow- 
ing we will mainly focus on the unitary regime where the 
equation of state takes a universal form. In the unpolar- 
ized phase the system is superfiuid at zero temperature 
and its equation of state is 
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where n = n-\ = n±, exhibiting the same power law den- 
sity dependence as in the noninteracting gas. The dimen- 
sionless parameter £ accounts for the role of interactions. 
Its value, evaluated with ab initio Monte Carlo simula- 
tions, is £ ~ 0.44 (see, e.g., yj). In the presence of po- 
larization we will assume that the system exhibits phase 
separation between an unpolarized superfiuid phase (s) 
governed by the equation of state J2J) and a fully spin po- 
larized phase (p) where interactions can be ignored and 
the equation of state is given by the noninteracting ex- 
pression 
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with n the density of the only species present in such 
a phase. This is the simplest assumption on resonance. 
The more complicated possibility of an intermediate non- 
superfluid phase where the two components coexist with 
different densities has been recently explored theoreti- 
cally 0,0 and experimentally 0. 
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In terms of the chemical potentials fi s and \i p the equi- 
librium between the pressures of the two phases can be 
usefully written in the form 



2fi s = (2£) 3/ % 
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It is useful to represent the equilibrium condition j2J in 
the [i— h phase diagram (see figure P). Here \i = + 
Hi)/ 2 is the average value of the two chemical potentials 
while h — — Hi)/2 is an effective magnetic field fixed 
by the difference between the two chemical potentials. 
From general thermodynamic considerations |5| we find 



(5) 



so that for h > 0, where jj, p = fi^, the phase separation 
between the superfluid and the spin polarized phase is 
described by the straight line 
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FIG. 1: Trajectories in the phase diagram for an unbalanced 
Fermi gas in a spin-independent potential (vertical dashed 
line) and for the dipole configuration discussed in the text 
(parabolic continouos line). The diagram exhibits a first order 
phase transition from the superfluid to a fully polarized phase 
along the diagonal dotted line. 
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The trapping introduces a position dependence in the 
chemical potentials of the two atomic spieces that can be 
evaluated according to the local density expressions 
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Result allows us to calculate the r-dependence of 
the chemical potential of both the superfluid and spin- 
polarized phases. The equilibrium condition |@J then be- 
comes an equation characterizing the r dependence of the 
surface separating the two phases. The constant values 
are determined by imposing the proper normalization 
on the the spin-up and spin-down densities. 

If the trapping potential is the same for the two species 
(Vj = V[) then the effective magnetic field h, which dif- 
fers from zero if ^ N±, is independent of position and, 
by varying r, we are consequently exploring the /i — h 
phase space at fixed h, as shown by the vertical line in 
Fig. 2] This is the case of all the configurations so far 
investigated in the literature. If instead, the two external 
potentials are different, the magnetic field h will be posi- 
tion dependent and, by varying r, we explore the phase 
diagram along the parabolic line shown in the same fig- 
ure. We will discuss in particular the case of a dipole 
displacement of two harmonic traps 
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with r± = ^/y 2 + z 2 . Furthermore we will assume that 
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the system to be globally unpolarized (iVj 
so that, for symmetry reasons, = /i" = hq 
magnetic field is simply given by the x-dependent ex- 
pression h — muj 2 xd. 

It is worth noticing that for any value of the displace- 
ment d the effective potential (Vj (r) + Vj. (r))/2 felt by the 
system in the superfluid phase has a minimum at x = 




FIG. 2: Sketch of a central slice of the dipolar configuration 
as given in Eq. <f 1 1 p . 



giving rise to a superfluid density profile symmetric in 
the x direction: 
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The potentials trapping the spin-up and spin-down com- 
ponents have instead minima at x = ±d respectively, 
thereby favouring the formation of a spin polarized con- 
figuration with density 
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The boundary separating the two density profiles is 
fixed by the condition Q of mechanical equilibrium and 
is characterized by the typical geometry of Fig. |2Jwhere, 
for sake of simplicity, we have shown a 2D cut. For exam- 
ple for x > we have fi p = /Uf and the boundary between 
the superfluid and the spin-up normal component is fixed 
by the condition 
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FIG. 3: Normalized induced spin dipole moment, D(d)/d, vs. 
displacement d of the trapping potential. The displacement 
is given in units of the Thomas-Fermi radius 

where we have introduced the radius R x = of 
the polarized component, A = lo±/lo x , a = (2£) 3 / 5 /(2 — 
(2£) 3/5 ) and 7 = 4((2£) 3/5 - l)/(2 - (2£) 3/5 ) 2 . We find 
that, for any value d < where R x = a l ^°(24N) 1 / 6 
is the Thomas-Fermi radius of the non-interacting cloud, 
with a x ° = y/h/mujx, there is equilibrium between the 
superfluid and the spin-polarized phases. Conversely, if 
d > R®, the superfluid is absent and the two spin polar- 
ized clouds are separated in space. 

Starting from results 10 and (|10|l it is possible to cal- 
culate the spin dipole moment Q as a function of the dis- 
placement d of the traps. The spin-up (-down) densities 
entering the integral are given by the sum of the su- 
perfluid and spin polarized components: n^i — n s + n^. 
Since the superfluid density is x-symmetric, it does not 
contribute to the integral and the spin dipole moment 
turns out to be reduced with respect to the value D = d 
predicted in the absence of interactions. In particular, 
when d <C R x the induced dipole moment takes the sim- 
ple power law behaviour 

2 10 / d \ 5/2 £ 9/S 
D{d) ^^\Rj) (2 ~(2C) 3 / 5 ) 5 / 2 ' (12) 

revealing explicitly that the magnetic response D(d)/d 
vanishes in the linear limit d — > 0. Figure |3] shows that, 
due to superfluidity, the induced dipole moment deviates 
significantly from the non-interacting value D = d at 
distances of the order of the size of the atomic cloud. 
At the same time the figure, as well as Ea. (|I2ll . shows 
that the induced dipole is always different from zero even 
for small displacements, revealing the absence of a true 
gap. These two fetaures characterize the behaviour of 
the trapped superfluid at unitarity. Actually, near the 
border of the cloud, the gap, being proportional to the 
Fermi energy, becomes smaller and smaller and even a 
tiny magnetic field can induce a finite, although small, 



magnetization. 

The polarization of the gas becomes more and more 
pronounced on the BCS side of the resonance. In fact in 
this case the gap decreases exponentially when we reach 
the border and, as a consequence, we can more easily 
polarize the medium by breaking Cooper pairs. For a 
trapped gas in the deep BCS regime we expect that the 
induced spin dipole moment approaches the value D = d 
at small displacement distances, which can be estimated 
using energy arguments as d ~ R X ^J A/ ep, where A is 
an average value of the superfluid gap. 

On the BEC side of the resonance, where molecules 
are formed at small density, we instead expect a very 
different behaviour, which can be described using 2-body 
physics. In fact in this case the gap remains finite even 
near the surface where it is fixed by the binding energy 
h 2 /ma 2 of the free molecules. The polarization will be 
zero until the value of the displacement is such that the 
molecular potential energy increase ~ muj 2 d 2 is of the 
order of the molecular binding energy. This corresponds 
to separation distances of the order of d/R° ~ l/k F a. 
In the deep BEC regime, where kpa <C 1, we therefore 
expect that the system will never be polarized unless the 
displacement is much larger than the size of the cloud [jj . 

In conclusion, the transition from the BCS to the BEC 
regime is predicted to reveal dramatic behaviour of the 
induced dipole moment which exhibits a transition from 
a quasi ideal regime where the gas is easily polarized and 
is only affected by superfluidity at small displacements 
(deep BCS), to a regime where the rigidity of molecules 
is strong and the system is not easily polarized by the 
separation of the trapping potentials (deep BEC). 

In Fig. |21 we have also shown the predictions for the 
induced spin dipole moment calculated using the phase 
diagram of the BCS approach which is known to pro- 
vide a semi-quantitive description of the equation of state 
at unitarity. In our approach the differences are due to 
the different value for the interaction parameter £ whose 
value is 0.58 instead of the correct value 0.44. The differ- 
ent curves reveal the sensitivity of the spin-dipole curve 
to the proper description of the equation of state. They 
are easily interpreted since the smaller is the value of £, 
the more favorable is the superfluid phase and, conse- 
quently, the smaller is the polarizability. The BCS ap- 
proach predicts that the phase diagram of figure ^ in- 
cludes a further, normal-mixed phase |(| which is simply 
a noninteracting gas of spin-up and spin-down compo- 
nents. At unitarity we then predict a transition from a 
superfluid to a mixed-normal phase and from a mixed- 
normal to a spin polarized phase as we move from the 
center to the border of the cloud. The inclusion of this 
additional phase, within BCS model, does not result in 
any visible effect in the physical quantity D(d). 

Separating the two trapping potentials and inducing 
a spin dipole moment is not the only way to polarize a 
Nf = Ni = N/2 system. Another procedure consists 
in modifying the trapping frequencies of the two spin 
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species, i.e, by choosing trapping potentials of the form 
V„(t) = \m{Ji^rl + < x x 2 ) , a =T, | , (13) 

with, in general ujj,* 7^ u;^,, for i = a;, _L. This gives rise to 
a relative compression of the two spin clouds. In the fol- 
lowing we consider the simplest case of isotropic trapping, 
for which the above procedure corresponds to inducing 
a spin monopole polarization where the radii of the two 
species take a different value. For a noninteracting gas we 
find R^-Rl = (a^°-a h l °)(24:N) 1 / 6 . In the presence of in- 
teractions the behaviour is quite different. In particular, 
even for large differences between the trapping frequen- 
cies the system remains fully superfluid and does not ex- 
hibit spin polarization unless the value of the interaction 
parameter £ is larger than 0.5. This behaviour is well un- 
derstood by a simple energetic argument. If the trapping 
frequency of one of the two spin species, say the spin- 
up component, is much smaller than the one of the other 
component, then the configuration with spin polarization 
would consist of two noninteracting, practically non over- 
lapping clouds with very different radii. The energy of 
this configuration would be EP/N = (3/4)fiw i (3A^) 1 / 3 , 
the contribution of the spin-up component being negligi- 
ble. Conversely, the superfluid, feeling a confining poten- 
tial with lo 2 s — {tJl + u>T)/2 ~ w?/2, has an energy equal 
to E s = This is smaller than EP since £ = 0.44. 

In other words, at least within the "two phase descrip- 
tion" , the system always prefers to form pairs and remain 
superfluid rather than giving rise to a phase sepration as 
happens in the dipole case. The effect of the asymmetric 
trapping is a change in the radius of the atomic cloud and 
in the collective oscillations as now the superfluid feels a 
different effective trapping frequency. Notice that this 
peculiar behaviour is typical of the unitary regime and 
of the BEC side of the resonance. On the BCS side one 
would expect a partial polarization of the system, since 
the energetic gain associated with the superfluid phase is 



smaller than at unitarity. 

The analysis of the dipole (and monopole) configura- 
tion has been here carried out at unitarity within the 
easiest assumption of phase separation between a totally 
polarized normal Fermi gas and a fermionic superfluid 
and with iVj = iVi . The generalization of our approach 
to iV-f 7^ N± configurations is straightforward. Further- 
more the inclusion of additional phases in the diagram of 
fig.l, resulting from the availability of more sophisticated 
microscopic theories, as well as of surface tension effects 
could be naturally accounted for in the calculation of the 
spin polarizability. 

Concerning the experimental feasibility of the mea- 
surement of the spin polarizability a conditio sine qua 
non is the possibility of producing spin-dependent trap- 
ping potentials. This in principle is feasible by profiting 
of the different polarization of the electronic spin in the 
two hyperfine states of the Fermi gas and working with 
magnetic gradients and/or polarized laser optical trap- 
ping (see [9j and reference therein) . However at unitar- 
ity - where the most interesting features concerning the 
BCS-BEC crossover take place and the superfluid phase 
is more easily achieved - the electronic spin polarization 
of two hyperfine states can be very similar, due to the 
large value of the magnetic field usually required to reach 
Feshbach resonances, making the experimental measure- 
ment of the spin polarizability a difficult task. This is 
indeed the case for 6 Li, while the situation looks more 
promising for 40 K, due to the smaller value of the mag- 
netic field at resonance. Let us finally remark that the 
formalism developed in this work can be easily general- 
ized to Fermi mixtures of different atomic species (e.g., 6 Li 
and 40 K). In this case the achievement of an independent 
tuning of the the two trapping potentials would be much 
easier. 
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